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SUMMARY 


A circular dielectric waveguide consisting 


of an isotropic core covered with a thin anisotropic 


sheet is considered. The sheet is represented as a jump 


immittance and Maxwell’s equations are 


applied. Solution of the boundary value problem yields the characteristic equation, or dtspers.on 
relation, whtch is then solved numeric*. The results are verified for the step-index fiber and area- 
ic wav.gu.des. Finally, examples are included to investigate the effects of the amsotrop.c 


lar, metallic 
sheet. 


vii 




1. INTRODUCTION 


1.1 Background 

Tie solution to the problem of were propagation in an isotropic, circularly symmetric optical 
waveguide is covered in many texts. In parfcular, each of references [lH 3] P-ides an analysis of 
the step-index opt.cal fiber. The step-index fiber is treated as a boundary-value problem wh,* ts 
solved by use of electromagnetic held theory and the separation of variables technique. The result 
is the characteristic equ.tton of the guide and its solutions yield the propagation constants of the 

modes existing in the waveguide. 

Recently, efforts have been made to And modal solut.ons for similar structures composed of 
anisotropic materials for apphcations in polarisation maintaining optical fibers. Related structures 
composed of metallic waveguides coated with anisotropic layers have also been analysed for radar 
eross section reduction purposes. In particular, Kapany and Burke [4) describe a circularly cyhn- 
drical waveguide with an isotropic cladding and an anisotropic core whose optical axis coincide, 
with the axis of the cylinder. Tonning (5) reformulates the problem using a matrix formalism, and 
then extends the work of Kapany and Burke by examining a waveguide with a uruaxtal core as well 
as cladding [6]. Another alternate approach is the variational analysis provided by Lind.ll and 
Oksanen [7] for the case of transverse anisotropy. Using this method, they compare then results 

the exact eigenvalue equation fo, a metallic, circular guide coated with a layer characterised 
by an isotropic surface impedance [8], For an anisotropic surface tmpedance, solutions have been 
obtained for propagation in the circular metallic guide |9[ . Chou and Lee analyse the propagatton 
and at.ennat.on character, sties of a multfiayered coated wavegu.de and then examine the radar 
eross section of the wavegu.de a. near axial incidence [10], Scattering from a multilayered cyhn- 
drical structure was studied by Gragba and Uslenghi [11] with each layer being represented by an 
anisotropic jump impedance. Thus there are a variety of methods available to analyse anisotrop.c 

structures. 


1 
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1.1 Outline of Research 

Using the jump immittance model, modal solutions will be obtained for a structure consisting 
of an isotropic, dielectric core covered with a thin anisotropic sheet surrounded by free space. 
The anisotropic sheet is modeled as a sheet of zero thickness with a jump admittance. The jump 
immittance model is convenient to implement and describes the electromagnetic behavior of a thin 
layer of penetrable material [12]. It is, therefore, possible to impose the appropriate boundary 
conditions and complete a field analysis using Maxwell’s equations. The derivation follows that of 
the step-index optical fiber boundary value problem. The only difference is in the application of 
the boundary conditions as dictated by the jump immitance condition. The resulting dispersion 
relation, or eigenvalue equation, given in Section 2.3 describes the mode structure of the guide. 
Also derived are the cutoff conditions of the various modes (Section 2.4). 

Since the jump admittance of the anisotropic sheet is generally complex in value, the dispersion 
relation and cutoff conditions are functions of complex variables. The solutions of these equations 
are the propagation and attenuation constants for the various modes that exist in the guide. As 
shown in Chapter 2, for limiting cases of the jump admittance, the dispersion relation and cut- 
off conditions approach the solutions for the isotropic step-index fiber and circular, metallic guide. 
To find the roots of complex- valued equations, the Newton-Raphson iterative procedure was em- 
ployed [13]. The FORTRAN program constructed to apply the Newton-Raphson technique to the 
dispersion relation is given in Appendix A. 

Chapter 3 provides the results of the numerical solution of the dispersion relation. First, the 
solutions are tested against known results for the step-index fiber and circular metallic guide. Next, 
cutoff values for the propagation constant are computed for the case of a lossless sheet. Finally, 
examples are included to illustrate how propagation and attenuation of the modes of the waveguide 
are affected by variations in the parameters of the anisotropic sheet. 
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1.3 Hpometry o f Waveguide 

The cross section of the circularly symmetric waveguide is shown in Figure 1. A dielectric core 
of radius p = a and index of refraction n = ^ is coated by an anisotropic layer. 



/O = CL 


Figure 1. Waveguide cross section. 


The layer is modeled with a jump admittance given by [H]- 


( 1 ) 


_ _ v , 1J°1 0 

110 “ y ° l o t?02/ u 


Yo = y/^o/fM) 


The coordinates of the sheet and rod are shown in Figure 2. 



Figure 2. Orientation of coordinate axes. 
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A coordinate transformation [14] is performed represented by: 


(2a) 

where 

( 26 ) 


*7 = R» 7 o rT - 


/ cos a 
y sin a 


— sin a 
cos a 


This transforms the arbitrary orientation of the sheet admittance into the p, 4>, z coordinates of the 
rod. The result is: 



with 

(4a) m = *701 cos 2 a -j- 7702 sin 2 a 

( 4 b) V2 = 7701 sin 2 0 + 7702 cos 2 a 


( 4 c) 


773 = sin a cos 0(7701 - 7702) * 




2. FIELD SOLUTION 


2.1 Field Components 

Maxwell’s curl equations for fields of angular frequency u are, in the absence of sources: 


(5a) 

(56) 


V x H = jure E 

V x E = H. 


We consider modes with longitudinal wave number f3 and attempt solutions of the form: 
(6) F = 

where F represents either the E or H field vector [6]. Note that 

8F 


( 7 ) 


dz 


= - 3(3 


Also, the periodicity of 2k radians in the variable requires that 

dF 


(8) 


d<p 


— jm\ m — 0, i 1, =t2, • ■ 


The curls of (5) axe expressed in the given cylindrical coordinates to produce a set of six scalar 
equations. With the use of (7) and (8), the transverse components E^.Ep.E^ , and H p may be 
expressed as functions of the longitudinal components E z and H z , as follows: 


(9a) 

- T"-> 

(96) 


(9c) 

1 — TTLuJ 6 dHr 

(9d) 

1 dE t m0 

Bt = 1; i - 3w ‘-d7 + T E - ) ' 


5 
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The transverse wave number k t is given by 

(10a) k 2 = k 2 -0 2 

_ ( = k 0 n, p < a 

(10^) — ( Wy/eojlo = ko, p > a 

For the region inside the core we must establish fields E z and H z which are finite at p = 0. 
Choosing 


( 11 ) 


k 2 = k 2 0 n 2 


we have for p < a 


( 12 ) 


F z oc J m (k t p)e 




Outside the core, an evanescent field is required, or 

k 2 t = k 2 -0 2 

(13) = < °‘ 


We have for p > a 
(14) 


F z « K^-upy”*-*', 


where K m is the modified Bessel function of the second type. This guarantees exponential decay 
as p - oo. Also, both k 2 and y 2 must be greater than zero for guidance. The bounds on (3 are then 


(15) 


k 0 n > (3 > ho . 
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The resulting modal components of all fields are: 
in the core, p < o: 


(16a) 

(16b) 

(16c) 

(16d) 

(16e) 

(16/) 


E z - AJ m (k t py m4, 30 ‘ 

H t = BJ m (k t p)e ,m4 ‘~ J0z 




a, = - t: B U i '^ ]e ” nt 


and outside the core, p > a: 


(17a) 

(17b) 

(17c) 

(17d) 

(17c) 

(17/) 


E, = 

lit = 

Ep = [ ^ CK . mM -^nK^)y^ 

JT, = l^2cK.b.0) - ^DK^(l,p))e’“*"’" 


2.2 A pplication »f Boundary Conditio ns 


oo as 


We have already satisfied the boundary conditions at n = 0 and ns p ■ 

, = a , the tangential electric field is continuous 


and (14). Across the interface at p 

(18a) 

(18b) 


(E l ) p =a. = (£*)*■=<*+ 
(E,p)p=a- ~ (^)p= a +' 


described by (12) 
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The tangential magnetic field is discontinuous and satisfies the jump condition [10]: 

(1®) P * (Up =a+ Hp=a+-) — ??(Etan)p=a » 

or, 

(20a) -J? z (a+) + H t {a _) = Y 0 [t^ E+ (a) + t) 3 E z (a)] 

(206) #*(*+) - ff*(a_) = y 0 [ij 3 tf*(a) + %£ x (a)]. 

Equations (18) and (20) will be applied to the field equations of (16) and (17). The result is a 
homogeneous system of four equations in the four unknowns .4, 2?,C,and D, 


( 21 ) 


M 


= 0 . 


For a nontrivial solution it is necessary that detM = 0. The resulting equation is known as the 
dispersion relation. 


2.3 Formulation of Dispersion Relation 


As shown in (21), the coefficients A, B , C,and jD,are chosen as the column vector and (18) and 
(20) are the rows. The right hand sides of (20) are evaluated at p = a_. The matrix M becomes: 





( 22 ) 


M = 


I J m (fcta) 

-YV m (M(r ?3 + m^) 


={2 jLM 

fct 


0 -KmM 0 \ 


Jm(fcto) 

-YoVl^ J lnM 


kfa 


-K m {l t a) 




\ + Y 0 (v^^mk t a) +yon.^i.(M) 




Taking the determinant of (22): 


(23) 




^j;(fc,a) 

' t 


J m {k t a) - Yo»h*ff J m( k * a ) 

0 

— K m (7( a ) 

%2-J m {k,a) + Yom lj ^ J ^{ k t a ) 

fc f a 



1 

iff Jm(M 



-y 0 Jm(fcta)(’?3 + T)i^) 


-Kfn^ta) 


^■4(M 

k, 


J m (k,a) 

-Yom^M 


r^-uw 

kfa 

jUfX 


-K m {y t a) 


= 0 . 


+ vo^ + ^ a ) J ^ kta) + Y ^-kT J - (kta) 


ZfrKmM 
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This last expression leads to the following result: 


n’Jfc l V X'mbfA-lelto + ?/*>) 

Ik? \k t J m {k t a)) ktlt k t Jm{kta) j'rtKmi'rta) 


k o 
7? 


klm y 

htK m {yta)) 




+ 


gVV. y, - „>) 

jfc t 7 ( 0 k t J m {k t a) jitKm{fta) 


u> 4(M A r ^(7tQ) 

+ Vtkt °^ 1 ktJm{kta) jliKmilt®) 


^ 0 n 2 fcg -Km(7t a ) \ 

ib ( k t J m (k t a) 7t htKmMJ 


Wy rji — 2— - 5 — J ?3 + T 72 ] 

+ T 0 W«MW # / 


7? a 


(24) 


a'MQy ^m(7to) + 2 ^tj 3 + »fc] = °- 

^7«**(WUW 


To develop a general ease, normalised variables [4] are introduced according to the following 


equations: 


(25a) 

bo = / 3 /fco 

(25i>) 

fcto = fct/ fc o 

(25c) 

0 

II 

0 

(25d) 

a 0 = afco 


The guidance condition becomes. 


( 26 ) 


n > bo > 1. 
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Equation (24) becomes: 

(JUM oooLV , JJJJwq) KM^° L( fr + 

£tMt VfctoJm(fc,oa 0 )/ JktoJ~(*«oao)7 t oJM7.oao) 




/ tf^.oao) N 

1 2 / mbp \ 

V7tO-Km(7tOOo)/ 

1 V ao / 

J^(fctoao) 

^(7to«o) 


JJM poo) ^(7toQo) / J ^( fct ° a °L 4- A ^( 7t ° a °L 

~ J>?lMr fc t0 Jm(^toao) 7toA r m(7toUo) V r fctoJm(Wlo) 7t0 m(7toao) 

. J^(fctoQo) rfjgM m _ 2-^-rj 3 + 172] 

T fctO^m(fctOCLo) \7t0 a 0/ 0tO a ° 


(27) 


+ 


^m(7toQo)_ [/ JgM r)i + 2^-% + rj 2 ] = 0. 
7tO-^m(7tO a o) UtO a o/ tO a O 


To write more compactly, let 


(28a) 


u =fcto a o 


(28b) 

(28c) 

(28d) 


q =r-y to ao 
uJ m {u) 

K M) 

m qK m {q ) 


Note that the dispersion relation is quadratic \xi J m - It is of the form 


( 29 ) 


CoJm + ^l^rn + C*2 — 0 ♦ 
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Collecting terms, 


( 30 ) 


JmCrM 1 - jm°oKm) 

+ + K m [/ir(l + VlV7 - vl) + «r] 

} 




+ ^i + 


<3-0 

jK m 

CLo 


/"m6 0 a 0 \ 2 mi> 0 ao 

I = — 1 TJI - 2 — r?3 + V 2 


V ? 2 i 

( mboao \ ' 

V “ 2 ) 


Vi + 2 


9 

m^oao 


i?3 + *72 




The step-index solution is contained in the real portion of (30) , while the anisotropic sheet adds 
additional complex terms. In this form, the dispersion relation is solved numerically; its roots are 

computed as explained in Appendix A. The solutions are [4]: 

_ -Ci + y/C? - 4C 0 C 2 

(31a) EH mn modes : Jm - 2 C 0 

_ -Ci - JC\ - 4C 0 Cz 

(316) HE ron modes : J m - 2 C 0 

Two cases are examined to test the jump admittance condition. The relative sheet admittances 
are complex: 

( g,,b t real, with: 

I g t — o lossless sheets, 

(32) t ?, = p< + ;6. < g t = b t = 0 no sheet present, 

( 1/6, = 0 perfectly conducting sheet. 

,. . r oisocst _ _ n obviously reduces to the equation for the 

Substituting the condition for no sheet, 17 , - 0, (3UJ odviousi.v 

step-index guide [1], 

(33) (<, + p.)-(”>W I (^ + p) =° 

Next, the metallic boundary of a perfect conductor is a particular case of an isotropic surface 

admittance („ = 0) in the limit as „ = * - oo. If (30) is multiplied by the factor the 

defining equation for the circular, metallic guide results: 


(34) ^(«)-U«) = 0 • 

At the limiting cases of no sheet and a perfectly conducting sheet, the derived dispersion relation 
matches the dispersion relations for the step-index fiber and the circular metallic guide. 
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2.4 Derivation of Cutoff Conditions 
2.4.1 Limits at Cutoff 

To find the cutoff conditions of the various modes, (30) is examined in the limit as 


(35a) 

(356) 

(35c) 


9 ~ * 0 
6 — * 1 

U — > aoytrMr — 1 


Under these conditions, K m takes the following form [4]: 


(зба) 

(збб) 

(36c) 


m =0; K 0 — f? 

ri n 2 
1 1,9 

m =l; + 


tti > 1 ; 


q 2 2 (m - 1) 

A detailed derivation of the cutoff conditions is presented in Appendix B. 


2.4.2 Case I. : m > 1 At first glance, the dispersion relation would appear to be a 

function of £ given the order of (36c). Upon combining the higher order terms and using the 
relation q 3 = ag(6g - 1), it is seen that ^ is the highest factor. Equation (30) is multiplied by the 
factor with the following result: 


(37) 





/ m 

3 m 

jv iMr 

— 



V a o 

1 f TTI 

( 

ao V a o 

m 


Go 


m - 1 


- 2jr) 3 H T - Mr (1 + VlV 2 - vl) - e r 


_ 1 (= - -2^1 - 2? - i (<=22)’,, + J=?« + „) = 0 

ao \ao m- 1/ it 2 a 0 V u 2 u ' 


A n um erical solution is computed from an equation of this form. The cutoff conditions are defined 
for modes in (31) with m > 1. 
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2.4.3 Case II. : m = 1 
by (36b). After performing a 


The m = 1 cutoff conditions are found using K x as defined 
similar simplifying combination as in Section 2.4.2, the dispersion 


relation is multiplied by Jm(u)j£r- The result is 


(38) 


Ji{u)[-jr)iHrCLoJ'i{u) + *M U )] - 0 


which implies the cutoff conditions: 

EH in inodes : Ji{u) = 0 


(39a) 

(39b) 


a 0 ^ 0 


HEi„modes : 


_ 1 

uJi(u) jv lMr a O 


2.4.4 r...m = 0 For modes of the order m = 0 ,.h« dispers,on relation factors 

easily. The resulting cutoff conditions from Appendix B are. 


(40a) 

(40b) 


TMon modes : Jo(u) = 0 

Jo( u ) _ ~ ^ 

TEo„ modes: - jWrMr 


Since the cutoff condttions (37), (39). and (40) are a specie! case of the dispersion relation, 
it is not surprising that they also produce the expected cutoff equations when examined for the 
limiting cases of no sheet and a perfectly conducting sheet. Note that the EH and HE modes of the 
step-index fiber become respectively the correspond^ order TM and TE mode, of the metalbc 
guide [10). For m < 1, the TM and EH mode cutoff remains unchanged regard],., of the lossless 
sheet parameters. The equation determining the cutoff frequencies becomes more complicated for 

higher order inodes. 

Given that the normalised propagation constant of (35) is real, the cutoff condttions are defined 
for propagating mode, which are lossless. Therefore, sheet parameters contamed in the precedmg 
cutoff equations me . .ken as pur, imaginary. Simple solution, for lossy sheets are difficult to obtain. 
However, such solutions are easily computed due to the availability of subroutines which calculate 

Bessel functions of complex arguments. 




3. NUMERICAL RESULTS AND CONCLUSION 


3.1 Lossless Sheet Parameters 

Due to the complexity of the dispersion relation, a numerical solution is attempted. Appendix A 
outlines the solution and provides the FORTRAN program constructed for the purpose of obtaining 
the zeros of the disperse relation. Lossless examples are examined first; the sheet admittance 
values are chosen to be pure imaginary. Marcuse [2] examines a step-index guide with «, = 1-0201 
and fj. r = 1-0 , Plotting normalized propagation constant against normalized radius. The results 
of Figure 3 agree with those shown in (2). The HE n mode is actually cutoff at zero, but the curve 
is nearly parallel to the a 0 axis at a 0 = 5 . The HE 21 mode is also plotted and it is not discernible 

from the TE 0 i and TM 0 i curve. 


o 



Figure 3. Normalized propagation constant 
vs. normalized radius for step-index fiber. 
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It, Figure 4, plots of the slep-.ndex guide ere provided with those of the corresponding mode, 
of the circular, met .Uic guide. Here, the horicontal axis has been changed to the cutoff number 
y 0oV 7^rrT. Remember, the TMo, cutoff value remains unchanged. In genereal, as Pn varies 
from -tO to -,oo (no sheet to a perfectly conducting sheet), the cutoff numbers will vary between 
the two curves shown. This variation is summarised in Table I„ the cutoff values are computed for 
some of the lower order modes. As the lossless sheet value, are decreased from sero, the step-index 
cutoff values increase towards the cutoff values of the circular, metallic guide. 



Figure 4. Normalized propagation constant vs. cutoff 
number, step-index fiber an d 

circular, metallic guide 





17 


TABLE I. CUTOFF VALUES 
FOR LOSSLESS SHEETS 
u = V At Cutoff 

T) ni, f?02» Q 

HEn 

HE 21 

TE 01 

-jo. ,-jo. , 

0 

0 

2.41 

2.41 

-j.2 ,-j.l , 

60 

.75 

2.58 

3.06 

-j.2 ,-j.l , 

45 

“86 1 

2.71 

3.14 

-j.2 ,-j.l , 

0 

1.02 

2.84 

3.26 

-j.3 ,-j.2 , 

0 

F24 

2.95 

3.41 

-i2 , -jl , 

0 

1.74 

3.05 

3.76 

-joo ,-joo , 0 

1.84 

3.05 

3.83 


3.2 T.ossv Sheet Parameters 


Suppose the sheet admittances are real or complex. Because loss terms are now present, the 
propagation constant will no longer be real. Instead, 


(41) 


1>0 = &0r + fo Os- 


The field equations will contain the term 


(42) 


e ko ' z e -jko ' z , 


where 6* is the normalized attenuation constant. Therefore, solutions with bo. < 0 will describe 
attenuated, propagating modes. 

The remaining graphs are plots of both the real and imaginary parts of the normalized propa- 
gation constant versus normalized radius. The normalized attenuation constants have been plotted 
up to the cutoff values of the propagation constant. Figure 5 represents the propagation curves for 
the TM 01 mode for sheet admittances that are real. The solid curve is the step-index fiber plot. 
Although the cutoff radius remained unchanged in the lossless examples, it has been increased for 
sheet parameters with small losses. For larger losses, the TM 0 i mode cutoff radius again becomes 

that of the step-index fiber. 





010 * 
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Figure 5. Normalized propagation constant versus normalized 
radius for TM 0 i mode with real sheet admittances. 


7701 ^02 Q 

0 o o 

.1 .1 0 

.2 .2 0 

1 1 o 
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The corresponding attenuation curves in F lg ure 6 reveal that a large attenuation near cutoff 
occurs for r*>, = -j.2 and becomes negligible near cutoff for no, = -H • The propagation constant 
curve for the no, = -jl case is not distinguishable from that of the step-index fiber. 

- 10' 5 



Figure 6. Normalized attenuation constant versus normalized 
radius for TM 0 i mode with real sheet admittances. 

tjoj ^02 a 

.1 .1 o 

2 2 0 

1 1 0 
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The same TMoi mode was examined for complex sheet parameters. While the propagation 
curves of Ftgure 7 are now nearly coincident, the attenuation near cutoff has been reduced signifi- 

cantly (Figure 8). 


o 



Figure 7. Normalized propagation constant versus normalized 
radius for TM 0 i mode with complex sheet admittances. 


r?oi 

7?02 

Q 

0 

0 

o 

■ l-j-2 


0 

-2-J.3 

-2-J.2 

0 

1-J2 

1-jl 

o 

(all four 

curves are 

nearly coincident) 




35.0 


4U.U 


T 


0.0 


5.0 10.0 15.0 20.0 25.0 30.0 


CL 


Figure 8. Normalized attenuation constant versus normalized 
radius for TM 0 i mode with complex sheet admittances. 


V02 ° 

,1-j.l o 

-2-j.2 0 
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The effect of real sheet admittances on the HE n mode is now examined 


Figure 9 shows the 


expected rise in cutoff radius. 



Figure 9. Normalized 
radius for HEn 


propagation constant versus normalized 
mode with real sheet admittances. 


T)01 V02 ° 

0 o o 

.1 .i o 

2 .2 0 

1 1 0 
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Note the slight change in scale of the attenuation curve of Figure 10. The lower values of Do, show 
less attenuation near cutoff than the TMoi mode. However, for the r, 0t = -Jl , the HE n mode has 
the greater attenuation. 



Figure 10. Normalized attenuation constant versus normalized 
radius for HE n mode with real sheet admittances. 
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Figure 11 shows agarn how the complex sheet parameters tend to move the propagation curves 
towards the metallic guide limit. 


O 



Figure 11. Normalized propagation constant versus normalized 
radius for HE U mode with complex sheet admittances. 
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Figure 12 shows that the complex admittances also reduce the attenuation of the HEn mode, but 
not as large a reduction as that which resulted for the TMqi mode. 


- 10' 5 



Figure 12. Normalized attenuation constant versus normalized 
radius for HEn mode with complex sheet admittances. 

Voi *702 Q 

•l-j.2 .1-j.l 0 

.2-j.3 -2-j .2 0 
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Finally, Figures 13 and 14 reveal the dependence of the propagation propagation and atten- 
uation curves on the sheet orientation. For the HE„ mode, increasing the parameter a has little 
effect on the cutoff radius, while slightly increasing the attenuation near cutoff. 



Figure 13. Normalized propagation constant versus normalized 
radius for HE n mode with complex sheet admittances 
and variation in the parameter a. 
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Figure 14. Normalized attenuation constant versus normalized 
radius for HE n mode with complex sheet admittances 
and variation in the parameter a. 
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The general effect of the anisotropic sheet admittance i. to raise the cutoff frequencies of the 
propagating modes. The cutoff frequncies vmy from those of the step-index Rber to those of the 
Circular, metallic guide as the sheet admittance vane, from aero to infinity. Use of lossy sheets 
also affected the cutoff of the lower order TM modes which were unchanged in the lossless case. 
Attenuated, propagating modes exist when lossy sheet admittances me present, although it was 
shown that the attenuation could be controlled by adjusting the sheet parameters and orientation. 

The dielectric waveguide analysed is a relatively simple structure. However, the method of 
analysis used is easily extended to mnltdayered structures. Additional boundary conditions will 
result, the jump immitance condition is applied at each sheet interface. The method of solution 
will remain the same; solutions must be found for a complex valued dispersion relation. 




appendices 


appendix a. dispersion relation root computation 

As stated previously, the dispersion relation, (30), is quadratic in J m . It is of the form 


(A 1) 


with, 


Co «7rn + C\J m + C 2 — 0 


(A 2a) 

(A2b) 

(A2c) 


C 0 = «rMr(l - ;m Q o^m) 

Cl = | — j^lMr a 0^m + £«n[Mr(l + T ^ lT ? 2 £f ^ 

1) 


ao 


(■mWj m _ 2 ^ va + m 


c >=^ +, v 


/ mi>o°o \ " 
/ 


m + 2 


mb o°o 


V3 + rft 


-{mb 0 ) 2 (^ + ^ 2 ) "° 


. 1 1 . f „_ XI tvn „ E-type modes are easily separated [3]: 

By applying the quadratic formula, solutions for H-type an . p 


(A3) 


F(a 0 , bo) = 2CoJm + Ci ± yjc\ - 4C 0 Cj 0. 


The HE modes are «he solut.ons taken from the positive radical and the EH modes are taken 
from the negative sign W . The method of solution is the Newton-Raph.on technique; i, is eapahie 
of determining complex root. |13] and thus accommodates both the iossless and loss, «ses.The 
dispersion relation shown in (A3) .« a funct.on of two tables. The Mowing FORTRAN program 
eolves this dispersion relation for selected parameters of the dielectric core and sheet admittances 
(cm.uo.,*,, -d «)• Initial jesses for the Newton-Raphson tecbmque are prov.ded from a 
knowledge of the cutoff conditions or from published data for the step-index fiber [2] . The calculated 
root i, then used as the next guess. A. one of the variables is incremented, this is continued to 

generate 
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APPENDIX A. (Continued) 

the dispersion and attenuation plot, over the desired range. Cutoff value, are obtained in a similar 
fashion, the cutoff equations are substituted in the function subroutine .ection of the program. 

****** * 

£ PROGRAM DISP * 

| \ 

1 in ill mm* hm I 

i r : rawaE «« * 

r RGES - INITIAL GUESS FOR BO * 

C FOUT : VALUE A 0F G d!IpERSI0N RELATION * 

l uS S S5 S d”(5evO, EB=1) : 

8 : 88 R al SlSi" 

c Sl 1 N2°N3 ANG - IWKeI PARAMETERS * 

8. **** 

C 

c 

T “» S,EKSE U1,E 1105 ' 

c BETWEEN SUBROUTINES. 

C COMMON /PARAM/ U.KTO Q.GTO N01 ,N02,ANG,EPSR,MUR.ITS2,NIN 

COMPLEX U.KTO, Q,GT0,N01,N02 
REAL ANG.EPSR.MUR.R 

INTEGER ITS2.NIN2 

C INTEGEr”mSTART~MSTOP~m”hSTARt”hSTOP ,H , ITS ,INC2 .AMIN 

REAL AO, AGES, REF, DEC 
COMPLEX BO.BGES.FOUT.F 
EXTERNAL NEWT.CBJNS 

C 3 FORMAT (IX, 12) 

C READ IN DATA AND INIT GUESS 

C OPEN (UNIT=7 9 , STATUS= ’ OLD ’ ) DIELECTRIC PARAMETERS 

C READ(79,7)EPSR,MUR 

7 FORMATClX ,2(F8 .4 ,2X) ) S HEET PARAMETERS 

° READ(79 ,9)N01 ,N02 , ANG 

9 FORMATClX ,5 (F8 .4, 2X)) MH RADIUS, MODE ORDER/TYPE 

C READ (79 , 1 1 ) AMIN , MSTART , MSTOP , HSTART .HSTOP 

11 FORMATClX, 5(14, 2X)) INIT . VALUES AND ITERATIONS 

C READ(79,13)BGES, AGES, ITS 

13 FORMAT (IX, 3 (F 10.5,21) ,13) 
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APPENDIX A. (Continued) 


C 

C 

c 


CL0SE(UNIT=79) 

OPEN OUTPUT FILE 

OPEN (UNIT=80) 

WRITE(80 , *) ( ’ EPSR 
WRITE(80,7)EPSR,MUR 
WRITE(80 , *) ( ’ NOl 

mstart hstop hstart 

WRITE (80 , 11 ) AMIN , HSTART ,MSTOP , HSTART , HSTO 
WRITE(80 ,*) ( ’ BO GUESS AO GUESS 

WRITE (80, 13) BGES, AGES, ITS 


MUR’) 

N02 


ANG’) 

HSTOP ’ ) 
ITERATIONS’) 


C 

C***CALC^ATC*ZEROS*OF*DISPERSION^R£LATION^^^^^^ |)(#i(i ^ i((i)ii)ii>i((i(tl(1<11(11(ll(1# ^ #<t< ,*#^ 

C *******WRlS(6 i)(’DO YOU WANT ITERATIONS PRINTED?’) 

WRITE (6 ',*)(’ REAL ROOTS: Y(l) N(2) COMPLEX :Y(-1) ,N (-2) ) 

WRITE(6 3 i)( ’INPUT RADIUS INCREMENT?’) 

WRITE (6,*) ( *DEC=IDEC+ . 1 ’ ) 

READ(6,3)IDEC 
REF=SQRT (EPSR*MUR) 

WRITE(80 , *) ( ’ MODE TYPE’) 

WRITE ( 80, *)(’H=’) ,H 
DO 1 H=HSTART, HSTOP 
DO 2 M=MSTART ,MSTOP 
BO=BGES 
AO=AGES 
R=FLOAT(AMIN) 

DEC=FLOAT(IDEC)*.l 
ITS 

WRITE(80,*)(’ ROOT FILE’) 

WRITE (80 , * ) ( ’ M= ’ ) ,M 
WRITE(80,*) (’ AO BO 
WRITE(80,39) 

39 FORMAT (3(/)) 


U 


FOUT’) 


DO LOOP SEARCHES FOR ROOTS WITHIH THE BOOHDS OF BO FOR 
PROPAGATING MODES 


14 

4 

2 

1 


DO 4 WHILE ( (REAL (BO) . GE. 1.) 

.AND. (REAL (BO) .LE. REF) 

.AND. (AO .GT. R)) 

CALL NEWT (BO , AO ,M,H ,FOUT) 
WRITE(6,*)A0,B0,U 
WRITE (80, 14) AO ,B0 ,U,FDUT 
FORMAT (IX, 2 (F10. 5) ,2X ,E12 . 5 ,4(F8 .4) ) 
AO=AO-DEC 
ENDDO 
CONTINUE 
CONTINUE 
CLOSE (UNIT=80) 

STOP 
END 



' COMPLEX FUNCTION F(B0,A0,M,H)^ 

"I — ... _ — , _ tmwm — — -» ~ ~ ^ * ** ~~ * 

_ . A Mtr 


?-j5^'iL5S'ssrf;uiS'i;s'S5insriiD trahsverse wave nos. 

c BETWEEN SUBROUTINES. 





ooooooooooooon 
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APPENDIX A. (Continued) 

COMMON /PARAM/ U KTO Q ,GT0 N01 ,N02,ANG,EPSR,MUR.ITS2,NIN2 

COMPLEX U,KT0 ,Q ,GT0 ,N01 ,N02 
REAL ANG.EPSR.MUR.R 

INTEGER ITS2.NIN2 


function f variables 

: ES (fs ARG * BESS FUNC.) 

KM - ABOVE RATIO FOR MOD BESS FUNC 

CO Cl C2 - COEFFICIENTS OF QUDRATIC FORM OF DISPERSION RELATION 
DISCR’ - QUADRATIC FORM DISCRIMINANT 


INTEGER H 
REAL AO, RANG 

COMPLEX BO *, JM KM , MB , CO , Cl , C2 , DISCR ,FK (5) ,F J (5) , C J 
PI=DC0NST( ’PI ’ ) 

N1~N01*C0S(RANG) **2+N02*SIN (RANG) **2 
N2~N01*SIN (RANG **2+NO2*C0S(RANG) **2 
N3* (NOl -N02)*C0S (RANG ) "'SIN (NANG) 

KTO=CSQRT (EPSR*MUR* ( 1 • , 0 . ) -B0**2) 
GTO=CSQRT(BO**2-(l . ,0.)) 

U=KTO*AO 

Q=GTO*AO 

CALL CBJNS(U,M+2 ,FJ) 

CALL CBESNK(GA,M+2,FK) 

IF (M .EQ. 0) THEN 

JM=(-1. ,0.)*FJ(M+2)/FJ(M+l)/U 
KM=-FK(M+2)/FK(M+l)/Q 

EL JM=(FJ(M)-FJ(M+2))/FJ(M+l)/U 

KM=-(FK(M)+FK(M+2))/FK(M+l)/Q 
ENDIF i ^ 

t -(MB*(1/U**2+1/Q**2))**2 

DISCR=CSQRT (Cl **2-4 . *C0*C2) 

IF (H .EQ. 1) THEN 
DISCR=-DISCR 

F=JM+ (Cl+DISCR) / (2 . *C0) 

RETURN 

END 


C 

c 

c 

Q 

s” ^*.::", ;”r;;y”So;wriiEsr-Ri?Bs5riTEM"m"rnCTSo^”r'* 

l cac 

C SUBROUTINE NEWT(BO ,A0 ,M ,H ,F0UT) 

MT ‘ im TRU,SVEESE U1VE ms ■ 

c BETWEEN SUBROUTINES. 

C COMMON /PARAM/ U.KT0.Q.GT0N01, N02.ANG.EPSR, MUR, ITS2.NIN2 

COMPLEX U.KTO.Q.GTO.NOl ,N02 
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APPENDIX A. (Continued) 

REAL ANG,EPSR,MUR,R 
INTEGER ITS2.NIN2 


SOBROmiNE KEJ^mliBLESj^^ op CIsreRSI oN RELATION OVER 
THE INTERVAL INC IN BO 


INTEGER M,H 

COMPLEX BO.FOUT.F.FP.INC 
INC=( . 10E-06,0 . ) 

IF (NIN2 .LT. 0) THEN 
INC*INC*(1.,1.) 

END IF 

l2 .RND U (RERLCBO)’ .LE E '(SORT(EPSR.mlR))) 

* F p=(‘(b6*inc ,io!« ,b”f (bo ,ro,h ,r) ) /i»c 

BO=BO-(F(BO,AO,M,H)/FP) 

FOUT=F (BO , AO ,M ,H) 

IF (ABS(NIN2) .EQ. 1) THEN 

WRITE(6 ,21) AO ,B0 ,U , Q ,FOUT 
21 F0RMAT(1X,9(F8.4,2X)) 

END IF 
1 = 1+1 
12 ENDDO 
RETURN 
END 





appendix b. derivation of cutoff conditions 
B.l Limits at Cutoff 

To find the cutoff conditions of the various modes, the dispersion relation 


" jmaoK-m) 


JJh 

ao 


f mh ° a A _2?^V3 + ri2 


( 51 ) 


+1C 2 m + 


{—) m 


is examined in the limit as 


(£2a) 


<7 — 0 


(B2b) 

(B2c) 


b - 1 


U — ♦ ao\Ar^f “ 1 


Under these conditions, K. m takes the following form [4]. 


(5 3a) 
(53b) 
(53c) 


m — 0 ; 

m = 1; 


m > 1; 


Ko 


K x 


q 2 l n ^f 


1 , !, 9 
— — -In- 
c 2 2 2 



B.2 Case I. : m > 1 

Equation (Bl) i. multiplied out and the factors of, am collected. Note that 


(54a) 

(54b) 



2 _ m 2 ™ 

*' Tn _ 9 4 (m - l)g 2 


+ 


1 

4(m — l) 2 
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APPENDIX B. (Continued) 

The dispersion relation will be of the fonn xtf + Y . Therefore, constant terms with no , 
dependence will be dropped. They will .ppronch aero ns the equation is multiplied by and , 
approaches eero. The surviving terms are 


q 

9 

- 2 ^ mb ° - m(M1 + " vl) + (t)] 


(B 5) 


- (— T )( fc 0 - *) (. 2 ) 


m 


2__ 2( =*£)»_ y™ [(^) ! n, + 2=ha + 1.1 
- 1 v u ao “ 


m 


u 2 


Since al(bl - 1) = 9 2 ,(B5) actually contains terms of highest order q \ , 


(B 6) 


JnffljVKr^ao 


- 2jT/3Mr - Mr(l + - *?2) “ M 

m - 1 / J 

m /m J1^3-i((-)V + 2 ! T'l3 + '») = 

"^aoVao m-lj “ 2 a 0 V u 2 “ 


u m a 


Multilplying (B6) by g 2 /m , the result is, 


JmiWrt^ao 

Jm \jVlH- 


(B 7) 


( HI !HL- 'j _ - Mr(l + V1V2 - 4) £t ] 

\ do TTl — 1 / 

1 (rn _ao_\?!H_i.((!HH£) 2 ni + 2^»)3 + %)=0 , 

- ao Uo rn- lj U 2 ao V a 2 “ 


which is equivalent to (37), the cutoff equation for the m > 1 modes 
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APPENDIX B. (Continued) 

B.3 Case IT, : m = 1 

For the m = 1 cutoff conditions, use K\ from (B3b) and 

K 2 _ I _ H + (inj) 2 . 

Kl ~ q A q 2 K 2 


(B 8 ) 


(Bl) is multiplied by , this allows 
i and constants with respect to <?• The surviving terms are 


removal of the terms containing <K„i' .terms np to the order 




• r ao Jm + l][Km - m26 ol = °- 


(B9) 

With appbcations of De l’Hospital’s Rule, this is equivalent to 

— jT)iiJ. r aoJi + 1 = 0 . 


(510) 

Equation (BIO) is the same as 
B.4 Case HI. : m =+Q 

Letting m = 0, the remaining terms of (Bl) are 


(39). the cutoff equation for the m - 1 modes. 


J 0 2 e r M r(l - JVi^o^o) 

J Mr ??2 \ 

-)-j7o| — JUiMtUO^-O + ^o[Mt(1 + VlV2 Vi) ao J 


(B 11) 


+K 2 + J — = 0 ■ 
&0 


This equation actually factors into the following two solutions 


(512a) 

( 5126 ) 


jc o(c f — nl) 

TEcase Jo = “ Mr (i _ J^K 0 

-1 /v . Vis 

TMcase Jo = + 


(5 13a) 


As <c„ becomes large at cutoff, these reduce to the quantities in (40) by inspection, 

TMon modes : *M U ) = 0 


(6136) 


J JM _ ^ ~ t h_ 

TE 0n modes : - jmao£rMr 


. . . . freauencies of the lossless TMoi modes are not 

From equation (B12b), it is apparent cutoff frequencies oi 


affected by sheet parameters. 
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